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Abstract 



This is the second paper of the authors in a series concerned with the development of a 
deterministic model for the transfer matrix of a MIMO system. Starting from the Maxwell 
equations, we have described in [I] the generic structure of such a deterministic transfer ma- 
trix. In the current paper we apply the results of [1] in order to study the (Shannon-Foschini) 
capacity behavior of a MIMO system as a function of the deterministic spread function of 
the environment, and the number of transmitting and receiving antennas. The antennas are 
assumed to fill in a given, fixed volume. Under some generic assumptions, we prove that the 
capacity grows much more slowly than linearly with the number of antennas. These results 
reinforce previous heuristic results obtained from statistical models of the transfer matrix, 
which also predict a sublinear behavior. 

Keywords: MIMO systems, Shannon-Foschini capacity, deterministic spread function, Fredholm 
determinants. 
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1 Introduction and the main result 

Multiple Input Multiple Output (MIMO) is a powerful technology for increasing data rates in 
wireless telecommunication. Experimental and theoretical studies show the increase of the ca- 
pacity (number of bits transmitted per second) when the number of the transmitting (TX) and 
receiving (RX) antennas also increases. Recall (see [ZIE]) that when the channel is unknown to 
the transmitter, the Shannon-Foschini capacity is given by 



where Mt is the number of TX antennas, Mji is the number of RX antennas, Imr is the Mr x 
identity matrix, Et is the average total energy transmitted by the TX antennas, Nq is the variance 
of the noise, H is the Mr x Mt channel transfer matrix which establishes the linear relationship 
between the signals at RX antennas and the signals at the TX antennas. 

Starting from the Maxwell equations, we have shown in [T] what is the generic structure of 
such a transfer matrix (see below (ll.2p ). In the present paper we apply the formula obtained in [T] 
and study the behavior of the MIMO capacity as a function of Mt, Mr and of the deterministic 
spread function of the environment. The antennas are assumed to fill in a given, fixed volume. 
According to [1], under certain conditions the transfer matrix can be well approximated by: 
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(1.1) 




(1.2) 



where §^ is the two dimensional sphere, s{Qii,flT) is a 6 x 6 matrix cahed the spread function 
which contains the scattering information of the environment, 



arinr) := {aff^nr) G C : ^ e {1, 2, . . . , 6}, m e {1, 2, . . . , Mt}} 

is a 6 X Mt matrix valued map which describes the radiation pattern of the transmitting system, 
while 

aRinii) :={a™^(r!fl) GC:me {l,2,...,Affl},je{l,2,...,6}} 

is a Mn X 6 matrix valued map which describes the receiving system. We assume that all a^" and 
a^'' are continuous functions of the angles. 

The index m accounts for the placement of the m'th transmitting/receiving antenna. Assuming 
that the transmitting/receiving antennas are placed in a finite volume Vt/r, the distance between 
them becomes smaller. Moreover, reasoning in terms of Riemann sums, we will assume that there 
exist two smooth kernels Al^^j^{fl, fl') such that: 



sup 



sup 



m—1 
, Ma 



m—1 



= 0{M^'), 

= 0{M^^), z,je {1,2,...,6}. 



(1.3) 



And here is our main result: 



Theorem 1.1. Assume that the spread function s{^r,^t) is the kernel of a Hilbert- Schmidt 
operator. We have the following situations: 
(i) . // Mji is fixed, then 

lim C{Mt, Mr) <oo; 



(ii) . // Mt is fixed, then 



^.^^ C{Mt.Mr) ^ Mt 



Mn^ca 



HMr) - ln(2) ' 

(iii). Let M = Mt and assume that Mr — aM for some constant a > 0. // the operator generated 
by the spread function of the environment has finite rank N < oo, then: 



C(M,aM) N 
lim sup — , . — < , . , . 
M^oo ln(M) - ln(2) 



(1.4) 



(iv). Let M — Mt and assume that Mr = aM for some constant a > 0. If the spread function is 
C°° in both angular variables, then for every e > we have: 



C(M,aM) „ 
lim = 0. 

M-J-oo M" 



(1.5) 



Remark 1. Let us go through some of the previous results obtained with probabilistic models 
for the channel transfer matrix. In the case when the distance between antennas is kept constant, 
some theoretical studies [H El E] conclude that the capacity grows linearly with the number of 
antennas. Still for probabilistic models, if the antennas are forced to occupy a fixed volume [71 18] 
then one has to consider correlations between them. This is done by introducing some ad-hoc 
correlation matrices depending on the interelement distances. In this case, they observe that the 
capacity either grows at most like a logarithm [5] , or even converges to some finite value [7] . In [5] 
it is also suggested that the capacity should tend to a limit as the number of antennas increases 
in a fixed volume. 
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Remark 2. In contrast with the probabihstic models where the correlations between antennas are 
introduced rather arbitrarily, in [I] we developped a deterministic ab-initio model for the channel 
transfer matrix (see (jl.2[) ') which implicitely takes into account these correlations, through the 
matrices ax and an which completely describe the radiation patterns of the transmitting and 
receiving arrays, while the spread function s{Q,n') describes the scattering environment. 

The mathematical technical assumptions in Theorem ll.ll are physically natural, thus our results 
confirm that the capacity of a system in a realistic environment grows more slowly than linearly. 
For example, in (i) we obtain that the capacity saturates when Mt grows and Mr is kept fixed; the 
physical explanation is that the spread function cannot convey enough transmit spatial diversity 
to the receiving side. Similarly, if Mt is kept fixed as in (ii), there is not enough transmit spatial 
diversity to start with and the capacity only increases as In(Mfl). When both Mt and Mr grow 
proportionally at the same time, then if there is not enough spatial diversity in the scattering 
environment as it happens in (in), we again only get a logarithmic growth. Finally, when both 
Mt and Mr grow and the spread function is varying smoothly, the growth is slower than any 
positive power of AI. 

Remark 3. Two important parameters which implicitely appear in (Hi) and (iv) are on one hand 
the value of the rank, and on the other hand the speed of oscillations of the spread function. Our 
proofs implicitely show that these factors are maybe more important in the capacity growth than 
the number of antennas. For the same distribution of antennas, the capacity should be larger if 
the environment contains a lot of scatterers and the spread function is very irregular. 

The structure of our paper is as follows: in Section [5] we express the capacity as a Fredholm 
determinant of an integral operator whose integral kernel depends on Mt and Mr in a way which 
is easier to deal with when these numbers grow. In Section |3] we prove that the capacity saturates 
as a function of Mt, while it can grow with Mr either logarithmically or power-like, but with 
arbitrarily small exponents. 

2 Shannon-Foschini capacity as a Fredholm determinant 

The main result of this section is contained in Proposition 12.41 but we need to start with a few 
technical lemmas. 

Lemma 2.1. Introduce the bounded linear maps At/r G i?([L^(S^)]^) defined by ("$ e [L^(S^)]^j; 




(2.1) 



Then At/r, are self-adjoint and non-negative operators. 



Proof Choose * = {V-^lLi ^ [L^{§^)f and $ = G [L'^{^^)f and compute: 





(2.2) 



and observe that At^) 



($,At^') = (At*,$). Moreover, 
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(2.3) 



The proof for Ar is similar. 



□ 
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Now let us consider B = HH* (the A//? x matrix appearing in the capacity formula (ll.ip ) 
and compute using (jl.2p : 



dflRdfl'j^ afl'(f2fl)6(i7ij, f2^)a^(r2^), where 



(2.4) 



We denote by B and S the integral operators in B{[L^{S'^)]^) given by the matrix valued kernels 
b{n, n') and s(0, f^') respectively. Note that the adjoint of S in B{[L'^{S'^)]^), denoted by S* , will 
have an integral kernel S*{i}, fl') = s*{n', fl). Thus from (I2.4p we have 



B = SAtS* 



(2.5) 



Lemma 2.2. Let K. := y/A^S* ArS^/A^ = {y/A^S y/A^)* {y/A^S ^/A^) > be a bounded 
positive operator in _B([L^(§^)]®). Then for any integer k > 2 we have 



B" = 



dQdn' aR{n)bk{^,n')a*ji{n'), 



where bk{^,^') is the integral kernel of the operator Bk = S\/ AtJC^ AtS* in _B([L^(S^)]^). 
Proof. Let us first show that the identity holds for fc = 2. We have: 



B' 



d^id^li^ aij(57i 



dn2dn^ b{ni,n2)a*R{n2)aRin3)b{n3,ni) 0^(^(4) 



diidn' aR{n)b2{n,n')a*j^{n') 

where b2{^R, ^'r) is the integral kernel of the operator 

B2 = BArB = SAtS*ArSAtS* = Sy/A^ICy/A^S* 

which proves the case k — 2. For k > 2 we have 
^k+i ^ ^k^ 



(2.6) 



S2xS2 



dflidfli aR{fli 



S2xS2 



dn2dn3 bkini,n2)a*R{n2)aR{n3)b{n3,^i) a*R{ni) 



S2xS2 



diidn' aR{n)bk+i{i^,^')a*R{n') 



(2.7) 



where 6fc(f2, fl') is the integral kernel of the operator Bk, and bk+i{^, is the integral kernel of 

Bk+l = BkARB = {Sy/A^IC''-^y/A^S*)AR{SATS*) = Sy/J^IC'' y/J^S* . 

The proof is over. □ 



Lemma 2.3. Let z belong to the intersection of the resolvent sets of B and K , i.e. z G p(B)r\p{IC). 
Then K, is trace class and we have the identity: 

TrcM, {{z - B)-^B] = Tr[i2(s2)]6 {{z - JCy^JC] . (2.8) 
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Proof. The operator K. is a product of two Hilbert-Schmidt operators, thus it is trace class. Since 
both sides of (|2.8p are analytic on p{B) D p{IC), it is enough to prove the equality for \z\ > 
niax{||i?||, ||A^||}. Using the power series expansion we can write: 

TrcM, ((z - B)-'B) = ^TrcM«E". 

n>l 

Then, using Lemma 12.21 and the trace cyclicity we get: 
Yl ^TrcM^B" = ^ ^Tr[i2(s2)]6{6„Afl} 

n>l n>l 

= E ^Tr[i2(s2)]45v/I^/C"-iv^5M^} = ^ ^Tr[i2(s.)]e {/C"} 

n>l n>l 

= Tr[i2(s2)]6{(z-/C)-i/C}, (2.9) 

□ 



Proposition 2.4. The capacity of our system can be written as: 

C(M.,M,).^Tr[,2(,2,eLn(l + ^/c), 

where JC = \/ AtS* AfjS\/ At , while Ln denotes the principal branch of the natural logarithm. 
Proof. The Shannon-Foschini capacity is equal to 

C(M.,M,) = ^Tr,..Ln(l + ^5 

Since the spectra of j^^B and jyf Mt ^ both included in some large enough closed interval 
of the type [0,L], we can find a positively oriented, simple and closed contour F which contains 
[0, L], while F C C \ (— oo, —NqMt/Et]. By the Dunford-Riesz functional calculus we may write: 

Ln f 1 + -^^b] =J- [ dz Ln(l + -^^z] {z - B)-^ 



Since the function g[z) := ^Ln [1 + j^ff^zj is also holomorphic on C \ (— oo, —NqMt/Et], we 
may write: 

Ln |1 + B] = — [ dz g(z)(z - B)-^z = — [ dz g(z)(z - B)-^B. 

Thus (using (gJl)): 

which can be integrated back and we obtain the result. □ 
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3 Proof of Theorem 11.11 

We start with an abstract technical lemma which will be used extensively during the proof. 

Lemma 3.1. Let F(x) = Ln(l + x) be defined on C\ (—00, —1]. Let Ti and T2 be two self-adjoint 
non-negative trace class operators defined on a separable Hilbert space %. Then F{Ti^2) o,re trace 
class and moreover, if AT :~ T2 — Ti we have: 

|Tr{F(T2)} - Tr{^^(ri)}| < ||AT||i, and (3.1) 
Tr{i^(r2)} < Tr{F(Ti)} + Tr{F'(Ti)Ar}. (3.2) 

Proof. If < i < 1 we define T{t) (1 - t)Ti + iTa- Clearly, T{t) > and T'{t) = AT € Bi{n) 
(the space of trace class operators). Let T a positively oriented simple contour contained in the 
analyticity domain of and surrounding the interval [0, max{(T(ri), i7(T2)}]. Then F completely 
contains the spectrum of T{t) for all t and we can write (in the sense of bounded operators): 

Fim) = ^ [ dz F{z){z - T{t))-^ = I - T{t))-^ - l/z}. (3.3) 

The second formula holds true because F{z)/ z is still analytic inside the domain of integration (see 
also the argument used in Proposition 12. 4p . But (z — T{t))^^ — l/z is now a trace class operator 
and it follows that F{T{t)) e Bi{H). 

Define the function : [0, 1] — > R given by: 

m Tr{F(T(t))} [ F(z)Tr{(z - Tit))-' l/z}dz. (3.4) 



27ri _ 

Let us first compute (j)'{t). Using the cyclicity of the trace we have: 

cb'it) = ^ ^ ^^(z)Tr{(z - Tit))-'T'it)iz Tit))-'}dz 

= 2^ / ^(^)Tr{(^ - Tit))-^T'it)}dz. (3.5) 

Denote by {\fj{t))}j>i the orthonormal basis consisting of eigenvectors of T{t), corresponding to 
the non-negative eigenvalues {Ej{t)}j>i counting multiplicities and arranged in decreasing order. 
Then using T'{t) AT we obtain: 

^'it) = Y,{f,{t),iAT)f,it))^ [ F{z){z-E,{t))-^dz 

i>l -^^ 

= Y.m)^ {AT)f,{t))F'{E,{t)) Tr{F'(r(0)Ar}. (3.6) 
j>i 

Thus (j)'{t) = Tr{(l + T{t))-'AT}, which means that |(/)'(t)| < ||AT||i and this proves ([511 . 
We now prove p.2p . Note the identity: 

0(1) = <^(0) + 0'(o) + / dt{\-t)4)"{t). 



Clearly, 0'(O) = Tr{T'(Ti)Ar}; thus the only remaining thing is to show that (j)"{t) < for all t. 
By differentiating once again in p.Sp we have: 

nt) = ^ ^ F'(z)Tr{(z - T{t))-'AT{z - T{t))-'AT}dz. (3.7) 
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Using the eigenvalues and eigenvectors of T{t) we get: 

-^"W^^EEK/^-W'^^^W)!' / F'{z)iz-E,it))-\z-E,it))-'dz. (3.8) 

j>l k>l 

Define {a,k}j,k>i, wliere a,fc(i) := F"{Ej{t)) if Ejit) = E^it), and ajfc(t) = ""''Iftijl^^tf; if 
i?j(t) 7^ Ek[t). Because F is concave on (0,(X)), all a^-fe's are non-positive. By the residue calculus 
we have: 

'/'"W = EEl(/jW'^^/fcW)l'«J'=W^O' 0<^<l> (3-9) 

j>l k>l 

thus the proof of p.2p is over. □ 



Remark. Assume that P = P* = is an orthogonal projection, and denote by Q = 1 — P. 
Define Ti = PT2P + QT2Q to be the 'diagonal' part of T2 with respect to the decomposition 
P + g = 1. Then T2 - Ti is off-diagonal and Tr{P'(ri)AT} = 0. Then jS^) implies 

Tr{P(T2)} < Tr{P(ri)} = Tr{P(PT2P)} + Tr{P(gr2Q)} 

which is a variant of Berezin's inequality. 

3.1 Proof of (i) 

The operators ^t/a (see (|2.ip ') have each a 6 x 6 matrix valued integral kernel Arp/p^{VL, 57') with 
the following structure: 

Mt 



m— 1 

Mr 

A)j(J7,l]')= E«™(^K"("')' »,je{l,2,...,6}. (3.10) 



A consequence of (jl.3p is the following estimate (in the sense of bounded operators generated by 
the corresponding integral kernels): 



||At/Mt-^t|| =0(Myi), llV^T/AfT- V^ll (3.11) 

Note that in general we cannot say more about the speed of convergence of the square root. 

Remember from Lemma that JC = ^/A^S* AnSy/A^. Identify EtK,/{NqMt) with Ti, and 
Et^/AtS* ArS^/At /No with r2. Since S is assumed to be Hilbert-Schmidt, we have that T2 — Ti 
is trace class with a trace norm which goes to zero with Mt- Thus p.ip implies: 



1 f Et 

m™ - ' ln-(2) ^'l^^f^^)!"^'' 1,^ + iVo 



\im C{Mt,Mr) = ^--^Tr[i2(s2)]6ln ( 1 + -L^s*ArS^] (3.12) 



and Theorcm ll.il (i) is proved. 

3.2 Proof of (ii) 

Let us introduce the operator 
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Another consequence of (jl.3[) and p.lOp is (again as bounded operators): 

sup \\Ar- MrArW <oo. (3.13) 

Ma>l 

Identify EtK:/{NoMt) with Ti, and MrVr with T2. Then \\T2 - Ti\\i is uniformly bounded in 
Mr, thus the estimate p.ip imphes: 



sup 

M„>1 



C(Mt, Mfl) - ^Tr[i2(s2)]6ln (1 + MrVr) 



< 00. (3.14) 



But now Vr is a rank A/j^, non-negative operator, and it is non-zero only on the range of Ay. 
Assume that Vr has exactly d positive eigenvalues, denoted by {Aj}^^]^, including multiplicities. 
Then 

1 1 

lim — — Tr[i2(s2)]6ln (1 MrVr) = lim V In (1 + AfflA^) = d<MT, 

Mn^oo \n(AlR) ^ ^ " Mr-^oo in(MR) 

and the proof of Theorem ll.il (ii) is over. 
3.3 Proof of (iii) 

Let us now examine the situation in which the operator generated by the spread function has 
finite rank. This would be the case if we have N isolated scatterers in the environment. Then the 
operator S whose kernel is given by the spread function can be written in the form 

N 

''J''\fj)(9k\, (3.15) 

where Cjk's are complex numbers, while {\fj)}jLi {|<7j")}|Li ^-^e (not necessarily unit) vectors 
in [L^(S^)]^. Remember that Mt = M and Mr — aM, with a > a constant. Let us introduce 
the operator 

Dm ■■= ^ y/Ar/MS* {Ar/M)S^J At /M , (3.16) 

thus: 

C(M, aM) = ^Tr[i2(s2)]6ln (1 + aMDu) ■ (3.17) 
The operator Dm can be written as: 

N I N N \ 



Dm = ^Y. {Ar/M) J2 CrnkU ) WAT/Mgj){^AT/Mgk\ 



j,k—l m— 1 

If we denote by: 

N N 



d,k := ( E^'j/^' (^r/M) c™*/™ ) ' l^j) - W^T/Mg,), (3.18) 



m— 1 



then 

N 



Dm = Y djk\hj){hk\. (3.19) 



We begin with a lemma: 
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Lemma 3.2. The matrix d is non-negative, and we have the following limit: 

N N \ 



lim djk =■ djU = — ^ ( Cijfi.AR Cmkfr, 



Proof. The convergence is implied by (I3.1ip , so we only need to prove non- negativity of d. Choose 

■1] = 

N N 

ih^d.^Lihi. — 

aMNo 



any {^/'jjjLi e and define J^Zi E.^i CzjV'j I/z) G [L^{S^)]'^. Then compute: 

d ^)c« = E E ^^^■fc'/'fc = ^«*) ^ 0, (3.20) 

j=i fc=i 

where we used the non- negativity of Afj, see Lemma [2. II □ 



Lemma 3.3. Let (pjk '■— {hj, hk) denote an N x N matrix constructed with the vectors introduced 
in (jS.lSp . Then (j) is non-negative, and we have the following limit: 

lim (t)jk ^- (l>jk ^ {gjjArgk}- 

Af — f oo 

Proof. The convergence is imphed by p.lSp , while the non- negativity of follows in the same way 
as for d, using the non- negativity of At proved in Lemma [2?T] We do not give further details. □ 



The following result expresses the capacity as a determinant of a finite rank matrix, uniformly 
in M. 

Lemma 3.4. We have the identity: 

C{M, aM) = j^ln det (l + aM Vd(j)V(tj . (3.21) 
Proof. It is enough to prove the equality 

Tr[i2(s2)]6 Ln (1 + zDm) = Trciv Ln (l + zVd(j)Vctj (3.22) 

for any z with \z\ sufficiently small. Then since both sides of p.22p are analytic functions in 
the half plane Re(z) > 0, the equality will also hold for z = aM. We will show that both 
sides of (|3.22p are given by the same power series around z = 0, which amounts to proving that 
Tr[^2(g2')j6 D^j = Tt£n {\/d(f)\fdY for ^-iiy P ^ 1- This is in fact a direct consequence of the 
identity (easily provable by induction) 

N 

Dli^ Y.mdY-%k\h,){hul p>l, 

in which one has to take the trace and use its cyclicity in order to move a \/d from left to the 
right. The proof is over. □ 



Lemma 3.5. The capacity grows at most like a logarithm: 



hmsup — ^TVT^ — < -^T^ < oo- (3.23) 
M^oo In(Af) - ln(2) 
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Proof. Seen as an operator on C^, \/d(j)y/d converges in operator norm to \/d (f) \/d when M grows. 
Due to regular perturbation theory, it follows that all N eigenvalues of \/d(j)\/d are uniformly 
bounded in M, say by a number A > 0. Thus 

C{M, aM) < ln(l + aM\) < {ln(M) + ln(l + aX)}, M > 1, (3.24) 

ln(2) lii(2) 

which ends the proof of Theorem ll.il (Hi). □ 
3.4 Proof of (iv) 

Remember that here we no longer demand S to have finite rank, but we assume that it has a 
smooth integral kernel. If Lb denotes the usual (non-negative) Laplace-Beltrami operator densely 
defined in i^(S^), then we denote by Lb ■= (B^^iLb the corresponding operator in [L^(S^)]^. The 
smoothness of s(fl, il') implies that for every natural number n, the operators 

are Hilbert-Schmidt operators. We know that Lb has purely discrete spectrum and the distribution 
of eigenvalues obeys Weyl's law. More precisely, for every E > define Pe and Pe to be the 
projectors corresponding to the eigenvalues of Lb and respectively Lb which are less or equal 
than E. Then it is well known [10] that: 

dim(FE) ^ ^ ^ 

E^oo E E^oo E 

The starting point of our proof are formulas (j3.16p and p.l7p . Introduce 

Ti^oMDm, T2 = aM^^AT/MS*PEiAB./M)PES^AT/M, (3.26) 

where T2 is obtained by inserting two projections Pe inside Dm- Then it is easy to see that there 
exists a constant Ci > independent of M and E such that 

||ri-r2||i<CiM||(l-PB)5||2, M>1. (3.27) 

11(1 - Pe)S\\2 = 11(1 - PE)Ls"Sn\\2 < -E-"||5„||2 

- r2||i < Ci M£;~"||5„||2, M > 1, E>0. (3.28) 

The crucial observation is that PeS has finite rank, equal to = dim(PE), and the method of 
Theorem ll.il (Hi) can be applied. The only problem is that we cannot be sure that the bound A 
in (j3.24l) can be chosen independent of N, but we will now show that in the worst case scenario A 
grows proportionally with N'^. 

Lemma 3.6. Let N ~ dim(PB). There exists another constant C2 independent of M and N such 
that: 

Tr{ln(l +T2)} < A^ln(l -l-CaMTV^). (3.29) 

Proof. We express Pe as X^jLi where the V-'j's are normalized eigenfunctions of Lb 

spanning the range of Pe- Then 

N 

i=i 

Comparing with p.l5p we see that fj — ipj, Cjk — 5jk and gj — S*^j. Looking at the definition of 
djk in p.lSp we see that \djk\ is bounded uniformly in j, k and M . Thus the norm of the matrix 



But 

for all > 0, thus 
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d can grow at most as N. The same conclusion holds for the matrix (jj defined in Lemma 13.31 It 
means that the norm of \fd(\)\fd (which is equal to its largest eigenvalue) grows at most as iV^. 
Thus we can choose some A = C-^N^ with C3 a constant independent of M and TV and use it in 
(P^ . The proof is over. □ 

Now using p.29p . (I3.28P and p. 251) in p.ip we obtain that for every n > 1 there exists a 
constant Kn > 1 independent of M and E such that: 

C(M, aM) < Kn {E\n{l + K^M E^) + ME^") . 

Now fix an < e < 1. Choose E = and introduce it in the above estimate. Since 

KnM E^ > 1, we have 

ln(l + KriM E^) = ln(l + R-Hr^ E-^) + ln(X„M E^) < ln(2) + ln(X„) + (1 + ^— ^) ln(M). 

n 

Thus we get another constant Kn > 1 such that uniformly in A/ > 1 and < e < 1 we have: 

C(M, aM) <Kn{l + Af ^ In(Af) + M"^ . 

But now we can choose rte to be the smallest natural number such that < e/2. Thus 
supj(j>j '^''^jlj-"^^'' < 00 and the proof of Theorem II. II is over. 

4 Conclusions 

In the case when the antennas occupy a given volume, our mathematical results will not constitute 
a big surprise for the engineers and researchers who have been involved in this type of MIMO 
studies and who have also predicted a sublinear behavior, even though they used ad-hoc statistical 
models for the transfer matrix. In this scenario, both our deterministic model and the stochastic 
ones seem to predict that the capacity growth can no longer be considered as linear if the number 
of antennas passes over a relatively low threshold. 

If the distance in between the antennas is maintained constant, the situation is rather different. 
All standard statistical models predict a linear increase in this case. But in a forthcoming paper 
we will apply our deterministic model in order to confirm the results of jllj which predicted a 
sublinear behavior even in this scenario. We will show that the sublinear growth begins from a 
not so large threshold value of the number of antennas, and that the reachness of the scattering 
environment is at least as important as the number of antennas. This shows that the discussion 
on the models is important. 
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